Theory of Ferromagnetism in Diluted Magnetic Semiconductor Quantum Wells 
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We present a mean field theory of ferromagnetism in diluted magnetic semiconductor quantum 
wells. When subband mixing due to exchange interactions between quantum well free carriers and 
magnetic impurities is neglected, analytic result can be obtained for the dependence of the critical 
temperature and the spontaneous magnetization on the distribution of magnetic impurities and 
the quantum well width. The validity of this approximate theory has been tested by comparing 
its predictions with those from numerical self-consistent field calculations. Interactions among free 
carriers, accounted for using the local-spin-density approximation, substantially enhance the critical 
temperature. We demonstrate that an external bias potential can tune the critical temperature 
through a wide range. 



The discovery of carrier-mediated ferromagnetisr 
in diluted magnetic semiconductors (DMS) has opened 
a broad and relatively unexplored frontier for both basic 
and applied research. The interplay between collective 
magnetic properties and semiconductor transport prop- 
erties in these systems presents a rich phenomenology 
and offers the prospect of new functionality in electronic 
devices. The possibilities for manipulating this interplay 
are especially varied when the free carrier system is quasi- 
two-dimensionaluTj because of the dependence of system 
properties on the subband wavefunction. In this article 
we address the dependence of the ferromagnetic critical 
temperature T c of two-dimensional DMS ferromagnets 
on the spatial distribution of magnetic ions, the subband 
wavefunction of the free carrier system, and on their in- 
terplay. Our analysis is based on a formulation of -the 
mean-field theory for carrier-induced ferromaffnetisrnB in 
a DMS which was developed in earlier workcl and is in- 
tended to be useful for any spatially inhomogeneous sys- 
tem. We find that T c is quite sensitive to the relative 
distributions of electrons and magnetic ions and can be 
altered in situ by applying a gate voltage. 

For definiteness we address the case of a (001) growth 
direction 2D hole gas in the Gai_ x Mn x As DMS sys- 
tem. Much of our analysis, however, applies equally well 
to other interfaces and to other DMS systems with cu- 
bic host semiconductors and, with one important caveat 
which we mention below, also to the case of a 2D electron 
gas mediated ferromagnetism. Our theory is based on an 
envelope function description of the valence band elec- 
trons and 9-spin representatioxL-of their kinetic-exchange 
interactionEl with d-electronsEJ on the 5 = 5/2 Mn ++ 
ions: 



H = U m + Hf + J P d Si ■ srfin - Ri), 
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the exchange interaction energy J p d ~ 0.15 eV nm 3 . In 
Eq.( |l|) TL m is the Hamiltonian of the magnetic ions, Tif 
is the four-band Luttingcr Hamiltonian for free carriers in 
the valence band, S is the magnetic ion spin and s is the 
electron-spin operator projected onto the j = 3/2 valence 
band manifold of the Luttingcr Hamiltonian. Wc assume 
here that the 2D free carrier density is sufficiently low 
that only the lowest energy heavy-hole subband is occu- 
pied, and that the quantum well of interest is sufficiently 
narrow that mixing between light-hole and heavvphole 
bands can be neglected. These simplifications leac£3 to 
a single parabolic band with the in-plane effective mass 
mt w O.llmo and the out-of-plane mass m* « 0.38mo. 
The two spin-states in this band have definite j z = ±3/2 
and definite s z = ±1/2; the projection of the transverse 
spin components onto this band is zero so that,thc kinetic 
exchange interaction has ran Ising character S3 

Our mean-field theorycl is derived in a spin-density- 
functional theory framework and leads to a set of phys- 
ically transparent coupled equations. In the absence of 
external fields, the mean polarization of a magnetic ion 
is given by 



(m)i = SBs(J pd S[ ni (Ri) - ni (Ri)]/2k B T) 
where B$(x) is the Brillouin function, 
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where i labels a free carrier, I labels a magnetic ion and 



~ 35 - 9053 -.*<!• (3) 

The electron spin-densities n a (r ) are determined by solv- 
ing the Schrodinger equation for electrons which experi- 
ence a spin-dependent potential due to kinetic exchange 
with the polarized magnetic ions: 

d 2 d 2 s i d 2 
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i>k,a{r) = e k ^k,a{r) ; (4) 
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In these equations v es (r) is the electrostatic potential, in- 
cluding band offset and ionized impurity contributions, 
v X c,cj(f) is an exchange-correlation potential on which we 
comment further below, and /(e) is the Fermi distribu- 
tion function. The spin-dependent kinetic-exchange po- 
tential, 



h p d(r) = J p d^S(r- Ri)(m) 



(6) 



is non-zero only in the ferromagnetic state. In the follow- 
ing we assume that the magnetic ions are randomly dis- 
tributed and dense on a scale set by the free carrier Fermi 
wavevector in the x — y plane, and that their density in 
the growth direction, c(z), can be precisely controlled.Q 
This allows us to take a continuum limit where the mag- 
netic ion polarization and the kinetic-exchange potential 
depends only on the growth direction coordinate and 



hpd(z) = J p d c(z) (m)(z). 



(7) 



In the following we discuss how the ferromagnetic tran- 
sition temperature T c depends on c(z). 

The T c calculation is greatly simplified when the spin- 
dependence of the exchange correlation potential is ne- 
glected and the quantum well width w is small enough 
to make subband mixing due to kinetic exchange inter- 
actions negligible, i.e., when Jc(z) <C ft 2 /m* z w 2 . We see 
later that the theory which results from these approxi- 
mations is normally quite accurate. With subband mix- 
ing neglected, the only effect of the mean-field kinetic- 
exchange interaction on the 2D carriers is to produce a 
rigid spin-splitting of the 2D bands, given by the first- 
order perturbation theory expression: 



Jpd / dz'(m)(z')c(z')\i/j(z')\ 2 



(8) 



where ip(z) is the growth direction envelope function of 
the lowest subband. For the case of a quantum well with 
infinite barriers, ip(z) — y/2/w cos(ir z/w). The electron 
spin-polarization can then be obtained by summing over 
occupied free-carrier states. It follows that, as long as 
both spin-f and spin-], bands are partly occupied, 
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Wz)\ 2 J pd I dz'(m)(z')c(z')\i>(z') 

(9) 



Eq. (|^) can be used to eliminate the free-carrier spin- 
polarization from Eq. (||) and to obtain a self-consistent 
equation for the function (m)(z). The equation de- 
pends only on the kinetic- exchange coupling constant 



Jpd, the free-carrier in-plane effective mass mjj, the sub- 
band wavefunction tp(z), and the magnetic ion distribu- 
tion function c(z). It is important to realize that the 
entire system behaves collectively. It has a single critical 
temperature which depends on the function c(z), rather 
than a z-dependent T c like that which appears in a naive 
adaptations of the bulk RKKY theory for T c to quasi-2D 
systems. 

An explicit expression for the critical temperature can 
be obtained by linearizing the self-consistent equation 
(^) . Expanding the Brillouin function to the first order in 
its argument gives the homogeneous linear integral equa- 
tion 



(m)(z) = 



where the kernel 



1 



k R T 



dz'K(z,z'){m)(z') 



(10) 



K(z,z') 



S(S 



l)J 2 P d m\\ 



12 



TV?? 



mz)\ 2 c(z')\i;(z')\ 2 . (11) 



Since the linearized equation is satisfied only at the crit- 
ical temperature, ksT c is equal to the largest eigenvalue 
of K{z, z'), the eigenvalue corresponding to eigenfunction 
f(z) oc \ip{z)\ 2 . It follows that 
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The critical temperature is proportional to J 2 d and m| , is 
independent of 2D carrier density, and for a system with 
uniformly distributed magnetic impurities is inversely 
proportional to the quantum well width, as illustrated 
by the solid line in Fig. |Tj. The accuracy of Eq. jl^ ) 
has been tested by numerically solving the set of self- 
consistent Eqs. (||)-(@) together with the Poisson equa- 
tion for v es (z) aad the local-spin-density-approximation 
(LSDA) formulali-3 for v XCtCr (z). If exchange and cor- 
relation effects are neglected (y XCtCr (z) = 0), the self- 
consistent numerical results represented by circles in 
Fig. |l| are identical to those obtained from Eq. (12) for 
quantum well widths w = 10, 15, and 20 ran. Interactions 
among the 2D carriers substantially enhance the critical 
temperature and also cause a strong dependence of T c 
on the 2D carrier density N, as illustrated by squares 
(N = 1 x 10 11 cm 2 ) and triangles (N = 0.5 x 10 11 cm 2 ) 
in Fig. |l|. Note that unlike the 3D case, where T c is 
propiactional to Fermi wavevector in the non-interacting 
limiCrca and is an increasing function of density even if 
the interactions are includedja the critical temperature 
for all three quantum wells of Fig. [j] is larger at lower 
density, M, = 0.5 x 10 11 cm 2 , where the interactions are 
stronger. I 
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FIG. 1. Ferromagnetic critical temperature of a quasi-2D 
DMS as a function of quantum well width w for fixed 3D mag- 
netic ion density c = lnm -3 ; inset: schematic diagram of the 
quantum well with the shaded portion of width d doped with 
magnetic ions. In this figure d — w. The non-interacting 
free carrier results were calculated both by solving the full 
self-consistent equations (circles) and by using the approxi- 
mate T c expression derived in the text (solid line) , confirming 
the accuracy of the latter approach. Squares and triangles 
represent the full numerical LSDA calculations for 2D car- 
rier densities N = 1 x 10 11 cm 2 and N = 0.5 x 10 11 cm 2 , 
respectively. 

Before discussing the implications of the expression 
( |l2| ) far T c we comment on its relation to the RKKY 
theoryu of carrier induced ferromagnetism in DMS's. As 
applied to bulk 3D electron systems the two approaches 
are fundamentally equivalent, although there are differ- 
ences in detail which can be important. The principle 
difference is that the RKKY theory treats the free-carrier 
magnetic-ion interaction perturbatively. As a result both 
approaches give the same value for the critical temper- 
ature where the magnetization vanishes; differences ap- 
pear at lower temperatures where the electron system is 
strongly or completely spin-polarized and a perturbative 
treatment of its interactions with the magnetic ions fails. 
Both approaches fail to account for the retarded charac- 
ter of the free-carrier mediated interaction between mag- 
netic ions in these low density systems, and for correlated 
quantum fluctuations in magnetic ion and free carrier 
subsystems which are important for some properties .EJj 
Both approaches are able to account for interaction ef- 
fects in the free carrier system, which increase the ten- 
dency toward ferromagnetism as illustrated in Fig. |l| and 
in the density-functional approach appear in the spin- 
dependence of the exchange correlation potential. It is 
in inhomogeneous cases, such as the reduced dimension 
situation considered here, that the density functional ap- 
proach has an advantage. Our approach provides a sim- 
ple description of the collective behavior of the quasi- 
2D ferromagnet and is able to handle geometric features 
which can be engineered to produce desired properties. 



For example we show in Fig. || a plot of the dependence of 
the critical temperature on the fraction d of the quantum 
well occupied by magnetic ions. When the integrated 
2D density of magnetic ions (c d) is fixed, the critical 
temperature increases as magnetic impurities are trans- 
ferred toward the center of the quantum well where the 
free-carrier wavefunction has its maximum. For fixed 3D 
density c the number of magnetic impurities decreases 
with decreasing d. This effect is stronger than the in- 
crease of the overlap between the carrier wavefunction 
and magnetic ion distribution resulting in the decrease 
of T c , as shown in the inset of Fig. |[ In both cases the 
analytic results of Eq. (jl^) are in qualitative agreement 
with the full numerical LSDA calculations. 
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FIG. 2. Dependence of T c on the central portion d of 
the quantum well (w = 20 nm) occupied by magnetic 
ions. The main graph is for fixed 2D magnetic ion density, 
c d = 20 nm -2 , while the inset is for the case of fixed 3D 
magnetic ion density, c = 1 nm -3 . 

One remarkable feature of quasi-2D ferromagnetic 
DMS's is the possibility of tuning T c through a wide range 
in situ, by the application of a gate voltage. In Fig. || we 
illustrate this for the case of a w = 10 nm quantum well 
with magnetic ions covering only a d — 3 nm portion near 
one edge. Even for such a narrow quantum well, the criti- 
cal temperature can be varied over an order of magnitude 
by applying a bias voltage which draws electrons into the 
magnetic ion region. 

In closing we remark that for several reasons DMS fer- 
romagnetism mediated by 2D electron systems is likely 
to occur only at inaccessibly low temperatures. One im- 
portant factor is the smaller value of the exchange inter- 
action parameter in the conduction band case. For ex- 
ample for the typical n- type, II- VI DMS quantum well 3 , 
(Zm_a;Cd x Se) m _/(MnSe)/jD the s — d exchange J sc i ~ 
4 x 10 -3 eV nm 3 . Our theory then gives T c ~ 10 niK for 
a 10 nm wide quantum well. Equally important, how- 
ever, is the weak magnetic anisotropy expected in the 
conduction band case, which will lead to soft spin-wave 



3 



collective modesS and substantial T c suppression. 
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FIG. 3. Dependence of T c on bias voltage applied across 
the w — 10 nm wide quantum well partially occupied by 
magnetic ions over a distance of 3 nm near one edge of the 
quantum well. Circles (squares) represent results of the full 
numerical self-consistent calculations without (with) the ex- 
change-correlation potential. 
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